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ABSTRACT 



SOME ASPECTS OF 
DEEPLY VIRTUAL COMPTON SCATTERING 



Elena Kuchina 
Old Dominion University, 2008 
Director: Dr. Ian Balitsky 

We consider different aspects of the virtual Compton amplitude in QCD on two 
examples: small-x physics accessible in the Regge regime and twist-3 approximation 
in the description of DVCS through the general parton distributions. Using this 
model, we give an estimate for the cross section of deeply virtual Compton scattering 
for the kinematics of CEBAF at Jefferson Lab. 
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CHAPTER I 
INTRODUCTION 

One of the inexhaustible studies in physics is the study of hadron structure . It is 
beheved that hadrons and nuclei are built from quarks and gluons, but it is still 
under investigation how those building blocks reveal themselves in nuclear reactions 
and how we could get information about those blocks from modern experiments. 

Theory offers Quantum Chromodynamics (QCD) as the gauge theory of strong 
interactions describing the dynamics between colored quarks and gluons. According 
to QCD, the interaction between the quarks becomes weak at very short distances. 
This phenomenon, known as asymptotic freedom, allows us to use perturbation the- 
ory to describe high energy strong interaction processes. At low energy, quarks and 
gluons are confined into colorless mesons and baryons which are the real world par- 
ticles. The QCD couphng constant grows at these energies plus the chiral symmetry 
of QCD is spontaneously broken by almost massless quarks. It brings us to the 
necessity of using non-perturbative methods at this regime. 

In search of ways describing the complexity of hadrons and hadronic systems 
one should perform precise low-energy experiments to gain information about the 
dynamics inside the hadron. However, one should also look at the specific high en- 
ergy kinematical regimes where the factorization theorems have been proven. Those 
experiments could reveal some new non-perturbative hadron structure information 
with the help of an accurate perturbative QCD description of the reaction dynam- 
ics [6-12]. 

The experiments performed with electromagnetic probes (Compton scattering) 
are the best ones. In Compton scattering, a real or virtual photon, emitted by 
a lepton, interacts with the nucleon with initial momentum p and as a result a 
real photon is emitted. Due to the new generation of electron accelerators and 
high precision, large acceptance detectors, very high precision Compton scattering 
experiments have become a reality. Despite the small cross sections of the process, 
it serves as a clean probe of hadron structure. 

It is convenient to examine this process in a parton model [13]. In this model, 
the nucleon consists of (a) three valence quarks, (b) quark-antiquark sea, and (c) 

The style specifications used in this thesis follow those of Physical Review D 
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gluons, the carriers of quark interactions. The sea and gluons is created and dis- 
sipated through time and is expected to have almost no influence on the process' 
outcome. The leading role in the nucleon belongs to the valence quarks. The typical 
time between interactions should be 1/A because A 200MeV is the only genuine 
scale in light quark QCD. The virtual photon, emitted by a lepton and carrying 
the momentum q {(^ — —Q^), has a lifetime 1/Q, which varies inversely with the 
momentum transferred by the photon. The momentum Q defines the virtuality of 
the photon. The higher virtuality, the shorter is a photon's hfe. At some point one 
may view the photon as being absorbed instantaneously by one of the quarks in the 
hadron. 

Suppose the quark that absorbs the photon has longitudinal momentum r. Upon 
absorbing the virtual photon, the struck quark becomes highly virtual with a lifetime 
r/Q^. Since this time is much shorter than the normal interaction time between 
quarks in the proton, the struck quark must reemit the photon before any interactions 
with the other quarks and gluons in the proton take place. 

Finally, since the transverse momentum of the absorbed photon is \ q \= Q the 
photon must be absorbed and re-emitted over a transverse coordinate region having 
I Ax I ^ 1/Q. The quark which absorbs the virtual photon is point- like (bare) down 
to a transverse size |A,t| ^ 1/Q. 

Thus, the scattering by the virtual photon takes place essentially instantaneously 
and over a very small, almost point-like, spatial region. Since the photon interacts 
only with a single quark and the transverse momenta r± of partons inside, the moving 
proton is assumed to be independent of photon virtuality Q^. It is expected that the 
amplitude of the process should be given in terms of the number density of quarks in 
the proton multiplied by the scattering amplitude of an individual quark. Charge and 
energy conservation brings us to the following sum rules for the parton distributions 
ni{x) inside the proton 



while the cross-section ai for the longitudinally polarized virtual photon is zero 




1 /-i 



(1) 



(2) 



and the cross-section (Tt for the transversely polarized photon is expressed in the 
impulse approximation as a sum of photon-quark cross-sections averaged over the 



distributions of quarks ni{x) and anti-quarks nj{x) in the proton: 

'^T^^e'^i^^x{ni{x)+nj{x)), (3) 
where Cj is the i-quark charge measured in the units of the electron charge e and 




FIG. 1: The diagram for deep inelastic scattering as an imaginary part of the virtual 
Compton amplitude. 

Most of the internal structure of the nucleon has been revealed during the last few- 
decades through Deep Inelastic Scattering (DIS). DIS is a type of inclusive scattering 
of high energy leptons on the nucleon in the Bjorken regime, which means that the 
photon virtuality is very large, xb — Q^/'^P • Q is finite, and the total center- mass 
energy of the photon- nucleon system s — {p + qY (one of the Mandelstam variables) 
is above the resonance region. 

The range of xb 

0<xb <1 

is given by the fact that the invariant mass of the unobserved final hadronic state is 
larger than the nucleon mass M 

{p + qf = + 2p- q + M'^ > 

and for the elastic scattering xb — i- 

In the process of deep inelastic scattering, the photon emitted by the lepton and 
absorbed by the proton creates an infinite number of possible final states. Summing 
up all those possible outcomes gives the intermediate state of elastic lepton-hadron 
scattering. This means that deep inelastic scattering could be described using the 
optical theorem by the imaginary part of forward Compton scattering. 
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Elastic scattering has the same initial and final momentum for both of the par- 
ticipating particles. The amplitude of this process would be dominated by short 
distance interactions in the Bjorken regime defined above. 

The leading contribution comes from so-called handbag diagram. The factoriza- 
tion theorem states that in the large-Q^ hmit the perturbatively calculable hard 
quark propagators completely factorize from parton distribution functions ni{x) 
{i — u,d,s, . . .), which describe non-perturbative long-distance information about 
hadronic structure. 

It was this particular process to which the perturbative quantum chromodynamics 
(QCD) was traditionally and successfully applied and where the approximate scahng 
behavior of the structure functions (independence of hadronic structure functions 
from the virtuality of the photon) was discovered. Unpolarized DIS experiments have 
mapped out the quark and gluon distributions in the nucleon, while polarized DIS 
experiments have shown that quarks carry only a small fraction of the nucleon spin. 
As a result, new investigations to understand the nucleon spin, both experimental 
and theoretical, became necessary. 

Yet, even more complex and intriguing information could be extracted from 
Deeply Virtual Compton Scattering (DVCS) [1,3]. As well as in DIS, the virtual 
photon carries large negative momentum q'^ = —Q^, but now the final momentum of 
the particles differs from its initial values. One should make sure that the momentum 
transfer t is as small as possible [1,2]. Again, the hard short-distance part factorizes 
from the non-perturbative long-distance part, general parton distributions (GPD), 
which now would be much more complex. 




FIG. 2: The diagram for virtual Compton scattering. 
In addition to usual variable x, the fraction of the nucleon momentum carried by 
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the photon absorbing quark, which runs from to 1 for the quark and from -1 to for 
antiquark, depends on two other variables. One of them is the momentum transfer 
t — {q — q'Y — {p' —pf which is an independent Mandelstam variable, that should be 
small and negative. The other one is the so-called skewedness ^, defined as a fraction 
of the average nucleon momentum carried by the overall momentum transfer in the 
process. These two extra degrees of freedom, t and ^, make the dynamics of deeply 
virtual Compton scattering rich and diverse [14-16]. 




FIG. 3: Parton picture for DVCS. Here P = {p + p')/2 is an average nucleon 
momenta, A — (p' — p) is the overall momentum transfer. 

Let us look at the handbag diagram of deeply virtual Compton scattering. Prom 
this figure it is easy to see that, depending on values of x and ^, the non-forward 
parton distribution functions can represent either the correlation between two quarks 
(when X > ^), two antiquarks {x < — ^), or between a quark and antiquark 
{—^<x<^). In the last case, GPDs behave hke a meson distribution and un- 
cover completely new information about nucleon structure which is inaccessible in 
DIS (which corresponds to the limit ^ — >^ 0). That is why DVCS is so compeUing 
and intriguing for physicists nowadays. 

There is only one disadvantage of measuring the virtual Compton scattering am- 
plitude, namely that the final photon can be emitted not only by the proton, the 
process we are most interested in, but also by an electron, which is referred to as the 
Bethe-Heitler (BH) process. This process is well known and it can be calculated ex- 
actly in quantum electrodynamics. The only thing which must be known is the clastic 
form factors of the nucleon. Unfortunately, light particles such as electrons radiate 
much more than the heavy proton. Therefore the BH process generally dominates 
the DVCS amplitude, especially at small t. 
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FIG. 4: a) the DVCS part of the amphtude; b), c) the Bethe-Heitler part. 

One way to minimize this problem is to find the kinematical regions where the 
BH process is suppressed (or at least is comparable with the DVCS amplitude). One 
of those cases will be discussed in Chapter 2. 

Another way is to exploit the interference between the DVCS and the calculable 
Bethe-Heitler process, by independently measuring both the real and imaginary parts 
of the amplitude. This approach requires the theoretical models which will allow 
one to calculate the interference term with high precision and better understanding. 
In Chapter 3, the leading (twist-2) approximation for the DVCS process will be 
discussed as well as the influence of the next to leading (twist-3) approximation on 
the flnal result. 

There are other promising non-forward hard exclusive processes such as longitu- 
dinal electroproduction of vector and pseudoscalar mesons at large Q^, but they are 
not discussed here. 
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CHAPTER II 
DVCS AT SMALL-X 

II.l INTRODUCTION 

In this chapter, the specific kinematical regime for the DVCS is considered when the 
energy of an incoming virtual photon Ee is very large p ■ g — > oo in comparison to its 
virtuality — oo, while the Bjorken variable xb = Q'^/2MEe is finite and small. 
It is the so called small-x region. To be specific, we calculate the DVCS amplitude 
in the region 

s > > -t > M^, (4) 

where s = (p + q)'^ — 2MEe is the Mandelstam variable corresponding to the center- 
of-mass energy squared of the photon-hadron system, M is the nucleon mass, t = {q~ 
q'Y is another independent Mandelstam variable equal to the momentum transfer. 
The mass of the lepton is neglected. 

The first study of the small-x DVCS was undertaken in Ref. [5]. The DVCS in this 
region is a semihard process in which quark ladders are dominated by gluon ladders 
well-known as the Balitsky-Fadin-Kuraev-Lipatov (BFKL) pomeron. It turns out 
that at large momentum transfer the coupling of the BFKL pomeron to the nucleon 
is essentially equal to the Dirac form factor of the nucleon Fi{t). Thus, the DVCS 
amplitude in this region ^ can be calculated without any model assumptions. Since 
there are only model predictions for the small-x DVCS in the current literature [19], 
even the approximate calculations of the cross section in QCD are very timely. The 
results obtained in this region can be used for the estimates of the amplitude at 
experimentally accessible energies where one or more conditions in Eq. Q are relaxed. 

To proceed further, at high energies it is convenient to use the Sudakov variables. 
Let us define 

P = P2 + —Pl, 

where {pi,P2) is the light-cone basis of the final particles momenta plane 

pI=pI = 0, 2pip2 = s ^ oo. (6) 



One advantage of these coordinates is their simple scaling properties when one takes 
the high energy limit. 
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All the other momenta are introduced in this basis as = akPi + l3kP2 + k± with 
Pik± = P2k± = by definition. For example: 

r = arPi + PrP2 + r_L. 

From = {p' — r)^ = and = {q' + r)^ = —Q^ one can estimate ar and f3r 
to be: 

^ -(Q^ + rl)/s. 
II.2 AMPLITUDE FACTORIZATION 

The amplitude of deeply virtual Compton scattering is determined by contracting 
the nucleon matrix element [20] of the T-product of two electromagnetic currents 
with the photon polarization vectors 

H = -^eMy;{q) J d'ze-^''-^p'\T{ji:,^,{z)j:„M}\p)- (§) 

It is known (see for example the review [22]) that in the leading order in perturba- 
tion theory the amplitude at high energy is purely imaginary up to the ^ corrections. 
At high orders in perturbation theory the amplitude will be purely imaginary in the 
leading logarithmic approximation (LLA) and one will restore the real part using the 
dispersion relations. 

Hence, we will first calculate the imaginary part 53 of the amplitude H 

V = -QH. (9) 

71 

The typical diagram for the DVCS amplitude in the lowest order in perturbation 
theory is shown in Fig. El (recall that the diagrams with gluon exchanges dominate 
at high energies). We are primarily looking for the imaginary part of the amplitude. 
At high energy the amplitude for the colorless particle scattering is described by the 
Feynman diagrams containing only two intermediate gluons with momenta k and 
r + k in the t-channel (Fig. E}. Simple estimations show that with good accuracy we 
can neglect the longitudinal momenta in their propagators: 

(r + A;)2 ~ (r + fc)i . (10) 
Let us consider the integral over gluon momentum d'^k = (Pk± i^^^Mk- 

^=11 TS^ I5FW^«* + -r,k), (11) 
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FIG. 5: The diagram for the virtual Compton scattering amphtude in the lowest 
order of the Leading Logarithmic Approximation (LLA). 

where $^^(A; + r, —k) and ($Ar)°^(— A; — r, k) are the upper and the lower blocks of the 
diagram with two gluon exchanges in Fig. IHl {a, b and ^, v are the color and Lorentz 
indices, respectively) corresponding to the impact-factor representation. 




FIG. 6: Factorization of LLA scattering amplitude into two impact factors. 

Due to our definition of the light-cone basis, all of the a-s in the upper block 
could be neglected as well as all of the /3-s in the lower block. At high energies, the 
metric tensor in the numerator of the Feynman-gauge gluon propagator reduces to 
9^^ f so the integral for the imaginary part factorizes into a product of 
two "impact factors" integrated with two-dimensional propagators 
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where 

Hk±,r^) = ^/2P2 j'^-^^Ku(k + r,-k) , (13) 

«fc=0 

and is the sum of squared charges of active flavors (u, d, s, and possibly c). 

II.3 PHOTON IMPACT FACTOR 

The photon impact factor is given by the two one-loop diagrams shown in Fig. [7| 
The first diagram yields: 

2 f d(3k f d'^p Tr[epp2{p + k)p2{p-r)t'{p- q)] 



sJ 2TTi J ?(27r)4 {p"^ + ie){{p - ry + ie) 

2m5{{q - pf)e{q^ - Po)2Tn5{{p + kf)e{po + k^). (15) 





FIG. 7: Photon impact factor. 

The second diagram differs only by trace and quark propagator: 
2 f dPk f d'^p Tr[epp2{p + k)e'{p + k + f — q)p2{p — ' 



X 



sJ 2m J 0271)^ {p"^ + ie){{p + k + r - qf + is) 

2m5{{q - pf)e{qo - Po)27iiS{ip + kf)9{po + ko). (16) 

In these formulae, the following notation was used: p = •y'^p^, S{{q — p)^) is a 
Dirac delta-function and 6{qo — po) is a Heaviside step function. The singularities in 
the physical region were replaced using the Cutkosky rule [23]. 

In a physical region the singularity is related to external 4-momenta. It is true for 
the normal thresholds. Anomalous thresholds do not show up as main singularities 
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in the physical region if the external particles are stable, which means that the 
4-momenta squared for each vertex is less than the smallest value of the normal 
threshold in the given channel. 

S. Coleman and R.E. Norton [24] found the following simple interpretation. In the 
physical region the main singularity of the diagram occurs only if all the vertices could 
be considered as some space-time points while the internal lines could be regarded 
as the trajectories of the real relativistic particles on mass shell. 

Due to the 5-function all the particles are on the mass shell 

P ^ H — 5 Pk ' 



as s aas 



and their energy is positive {9{q — p)o 9{q — p,p2) 9{d^) and 9{p + A;)c 
9{p + k,p2) 9{a^)). Here we introduced the notation a = 1 — a. 
The trace in the integral ()16|) contains fc-dependence: 

—^Tr[epp2{p + k)e'{p + k + r - q)p2{p - q\ = 



-(e,e')2(pl + i?xP±) 

-(e,rx)(e',^)4pi 
s 

+ (e,i?^)(e',px)2 

+ (e, ^)(e',pi)4((l - 2aa)pl + R^p^) 
s 

-(e,^)(6',^)8pi(pi-r^px) 
s s 

-(6,^)(6',px)4pi(l-2a) 
s 

-(e,^)ie',R^)Apl 

-(e,P±)(e',pi)4aa(l - 2a) 

-(e,p^)ie',^mi-2a)pl + 2R^p^) 
s 

+ {e,p±){e',pA_)8aa 

-(e,P±)(e',i?±)2(l-2a) 

-(e,j9i)(e',pi)8a^a^ 

+ (e,pi)(e', ^)4((1 - 2aa)pl + (1 - 2a)R±pi_) 
s 

-(e,pi)(e',p^)4aa(l - 2a) 

~{e,pi){e',R±)4:aa 

+ i^,R± + ar^){e',^)AR^p^ 
s 
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+(e,P±)(e', R± + oir±)Aa{l - 2a) 
+(e,pi)(e', it!_L + ari_)%a'^a 
+ (e, — )(e', -R_L + ar_L)8ap^ 

-(e,rx)(e',^)4(i?^r^ + arl) 

-(e,p^)(e', ^)4(i?^ + ar^)((l - 2a)i?^ - a(4p^ + (1 - 2Q;)r^)) 

— (e,pi)(e', — )8Q!(i?_L + ar_L)((l — 2a)p± + aR± — aa.r±) 
s 

-(e, ^)(e', ^)8(i?x + arx)(2p^ + R^- ar^)pl, (17) 



hidden in the definition of R±: R± = k± + ar±. 

This formula is rather long, but its beauty is in its generality. Here one has no 
restrictions on photon polarizations or virtualities. It simphfies a lot for every specific 
case. For example, when both photons are transverse (e,pi) = (e,P2) = = 
(e',p2) = 0, this integral reduces to: 



The integral for the first diagram differs only in overall sign and the definition of 
R± = —ar±. By introducing the Feynman parameter a' (with simultaneous shifting 

to p± ^ p± + a'R_L) and using formulae of Feynman parametrization, one is able 
to get the photon impact factor for the case of arbitrary photon polarization as the 
following: 



s 



s 




x(-2) (e,e')(pi + i?±p±) 

— {£,P±){^' ,r±)2aa{l — 2a) 

- (e,p± + R±)(e',p±) 

+ (e,p^)(e',i?^)(l-2a)2 . 



(18) 




(e, e')(aa(l - 2aa)(l - 2a'){q''^ - q^) + (1 - 2aa - 2a'a' + d>aaa'a')R\)) 

■(e,pi)(e',pi)8a^a^ 

(e, —){e' ,pi)2(aa{l-2aa){l-2a'){q''^ - q^) + 
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(l-2aa-2aV+4aaa'(2-a'))^i) 
+ {e,R±){e,pi)Aaa{l - 2a)a' 

-(e,pi)(e', y)2(aa(l - 2aa)(l - 2a'){q'^ - q^) - Aa^a^r\ + 
4aa(l - 2a)a'(i?_L, r^) + (1 + 2aa - 2a'a' - 4aaa'^)i?y 

-(e, —Me', —)4(aa(q'^ + rl)(aa(l - 2a')(q'^ - q^) + (1 - Aa' + 2a'^)Rl 
s s ^ 

(1 - 2a){R±,r±){aa{l - 4:a'a'){q'^ - q^) + (1 - 6a' + lOa'^ - Qa'^)Rl 
-(e, i?^)(e', ^)4aaa'((l - 2a)(g'2 + rl) + 4a'(i?^, r^)) 
+ (e,rx)(e', — )2faa(l - 2aa)(l - 2a'){q'^ - q^) + 

(1 — 2aa — 2a' a' + 8Q;aa'a')i?5_) 
-(e,pi)(e', i?_L)4aa(l - 2a)a' 
+ (e,pi)(e',r^)8aV 

+ (e, ^)(e', i?^)2(l - 2a)(aa{l-2a'f{q'^-q^) + (l-6a'+10a''-6a'3)i?l) 



-(e, — )(e', rx)4aa(aa(l - 2a')(g'^ - Q^) + (1 - 4a' + 2a'^)R 
s 

— (e, i?^) (e', i?j_)8aaa'a' 
-(e, i?^)(e',r_L)4aa(l - 2a)a' 

— R^ = + ar± —>■ R!^ = —dr± 



,/2\ d2 



n'^^(a'q" + a'q'). 
a'a' 



(19) 



(20) 



One can see from Eq. ()19p. that even in the most general case, the difference 
between the two diagrams of Fig. U\ is in the presence of k± in one and the absence 
{kj_ 0) in the other. It means one can write [25]: 



Hk^,r^) = l{k^,r^)-I{0,r^). 



(21) 



X 



For the particular case of two transverse photons one gets (cf. [26]): 

|(l-2aa)/?l(e,e')x 
+4aaa'[i?i(e, e') - 2(e, R)±{e', R)±\ 
-4aa(l - 2a)(r, e)±{R, e')±\, 



(22) 
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and for the longitudinal polarization 

1 



e\q) = ^{pi + xp2) (23) 



of the incoming photon the formula is also simple: 

^ ' ^ 2Q Jo 271 Jo 271 ^ ^ ^ 'i 
|(l-2aa)i?i(r, e')± 

+4aaa'[Rl{r, e')± - 2(r, R)±{e', R)±] (24) 
-4aa(l - 2a)Q2(i?, e')±|. 

Here (a, 6)_l denotes the (positive) scalar product of transverse components of 
vectors a and b. At large transverse momenta k'j_ ^ r'j_, the impact factor (j2H) 
reduces to: 

The impact factor for the proton, which describes the pomeron-nucleon coupling, 
cannot be calculated in perturbation theory. However, in the next section it is 
demonstrated that at high momenta k'j_ ^ this impact factor reduces to 



/iv(fcx,rx)''=''Vr"(t), (26) 
where Fl^^{t) is the sum of the proton and neutron Dirac form factors. 

II.4 NUCLEON IMPACT FACTOR 

In the lowest order in perturbation theory there is no difference between the diagrams 
for the nucleon impact factor shown in Fig. |Hland similar diagrams with two gluons 
replaced by two photons. One has to add the trivial numerical factor cp = ^ and to 
replace e ^ g. In this case the lower part of the diagram can be formally written as 
follows: 

$^(_A; _ r, k) J dze''^p'\T*{J,{z)MO)}\p), (27) 

where = u'^^u + d'-j^d. The T* means that we take into account only the T-product 
of the diagrams with pure gluon exchanges in the t-channel excluded. By definition, 
such diagrams contribute to subsequent ranks of the BFKL ladder rather than to the 
impact factor. This is also the reason why we have not included strange quarks in 
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FIG. 8: Nucleoli impact factor. 



the definition of electromagnetic current J . Since fc^ in our case is large and negative 
(-fc2 = k\'^ M^) we can expand the T-product of two currents near the light cone 
(see e.g. [27]) 



where again (...)* stands for the matrix element with pure gluon exchanges ex- 
cluded. Here [x, y\ denotes the gauge link connecting the points x and y 
{[x,y] = Vexp (ig Jq du{x — y)^A^{ux + (1 — u)y)j). The matrix element can 
be parametrized in terms of skewed parton distributions [28] as follows: 

{p',X'\q{z)[z,0]p\q{0)\p,X)U^, = u{p' , X')p\u{p, X) /We^(^-^)^X'(X, t) 

Jo 

+^«(p', X')p\fMP, A) j^We^(^--)^^W,«(X, t), 



{p',X'\m[0,z]p\q{z)\p,X)U^, = u{p\X')pMvA) / dXe'^''^^Vl{XA) (29) 



^-u{p\ X')prMP, A) f^dXe-^^'^^WliX, t), 



where V" (X, t) and W" (X, t) are the non-flip and spin-flip skewed parton distri- 
butions for the valence u quark (recall that one must take into account only va- 
lence quarks since diagrams with pure gluon exchanges are forbidden). Similarly, 
V^(X, t) and W^(X, t) refer to the valence d-quark distributions. At large energies 
u{p', X')piu{p, X) = s6xx', so 

{p', X'\m [0, z]p,q{z) - q{z) [z, 0]pig(0) |p, X)*,,^^ = (30) 



dX U 

I ^ 



-iXpz i(X—x)pz 



s6,yV',iX, t) + u{p\ A')^«(P, A) W,^(X, t) 
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After integration over z the lower block ()27p reduces to 



$Ar(-A; - r, k) 



dX 



(31) 



3s Jo 



(X - x)s + 2pi-k 



-Xs + 2pi- k 



-k^ - 2p ■ k{X -x)-te -k^ + 2p-kX -ie 



[SMV:iX,t) + V'AX,t)) +u{p\X')^^u{p,X){W:{X,t) + W'AX . 

The nucleon impact factor ()14|1 is the integral of the imaginary part of right hand 
side of Eq. (13 Ij) over energy: 

,2 



lN{k±,r±) 



1 dak^^ , I rl 
-— b<I>Af(-(aA: )V\->^L-r^,(XkPx^kx 

ZTT S 



(32) 



1 .1 .1 
- (iafc / dX 



s[X - x)b{k\ - OLks{X - x)) - sXb{k\ + (XksX) 



S,y{v:{X,t) + Vt{X,t)) +u{p',\')^^u{p,X){mx,t) + W',{X,t)) \ . 



And finally, 



lN{k^,r^) = y^'dx(^6,y{V:{X,t)+V',{X,t)) 



(33) 



2Ms 



u{p\ x')pifMp, A)(w:(x, t) + W^{X, t)) . 



Since valence quark distributions decrease at x — one can extend the lower limit 
of integration in the r.h.s. of Eq. (j32|) to and obtain: 



lNik±,r±) 



+ 



3 
1 



^'rfxkv(v:(x,t) + v,^(x,t)) 



2Ms 

Let us recall the sum rules [1,28 



uip',X'W±u{p,X){W:iX,t) + W',iX,t))]. (34) 



rfX(J-,^(X,t)-^|(X,t)) = Ff(t), 
dX{ICl{X,t)-ICl{X,t)) = Flit), 



(35) 



where J^^{X,t) and /C^(X, t) are the total (valence + sea) non-fiip and spin-fiip 
skewed quark distributions while J-'^{X,t) and /Cf(X, t) are the antiquark ones. In 
these equations, Fi{t) and Fl{t) stand for the g-quark components of the Dirac and 
Pauli form factors of the proton. Since the contribution of sea quarks drops from 
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the difference T'^ — T"^ one can rewrite Eqs. ()35p as the sum rules for valence quark 
distributions: 

C dXVl{X,t) = F?{t), f'dXWl{X,t) = Flit). (36) 

Jo Jo 

Substituting this estimate to Eq. (jH^ and using the isospin invariance, one can get 
the final result for the nucleon impact factor at large transverse momenta 

Mk^, rx) S,yFnt) + i^u{p\ X'WMp. m^'^it), (37) 

where Ff+"(t) = Ff{t) + Fi"(t) and ^^"(t) = F^{t) + F^{t). As usual, Ff^"^ and 
Ff*^"''* are the Dirac and Pauli form factors of the proton (neutron), respectively. With 
our accuracy they can be approximated by the dipole formulas: 



Fl + Fit/AM' =G\ =l/(l + |t|M)2, 

+ =G^^ = 2.79/(1 + |t|M)2, 

+ F^t/AW^ =G% =0, 

Fl^ + F^ =Gl, = -1.91/(1 + |t|M)2, 



(3^ 



which leads to 



Ff+"(t) 



1 l+0.88|t|/4M^ 



1 



irm^ 1+I*I/4M2 (i+|t|/a,)2' 



Fr\t)= ^0. (39) 

Here the notation at = 0.71 GeV^ is introduced. In further calculations we will use 
the estimate (jSH). It is obvious that the spin- flip term turned out to be negligible for 
our values of t. Besides, it vanishes at t = 0, which suggests that it is numerically 
small at all t. 

Our final estimate of the nucleon impact factor is: 

lN{k±,r^) 5AA'Fr"(t), (40) 

where Ff"*"" is given by the dipole formula ()39p . 

The dipole formula for the neutron form factor does not seem to work as well as 
the dipole formula for the proton form factor. As a measure of the uncertainty one 
can compare the results obtained from Eq. (jSHl) to those obtained using the model 
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from Ref. [29] (which was fit only to the proton form factor): 

Fnt) = l£dx{v:{x,t)+v'Ax,t)), 

V^{X,t) = 1.89X-°-^X=^-5(l + 6X)exp ^ 



V X2.8GeVV 

V'^{X,t) = 0-54X-°-^^^-^(l + 8X)exp(^-|^^A^). (41) 

The results for the DVCS cross section in this model are about 1.5 times bigger than 
the results obtained from the dipole formula (jH^ . 

In what follows the factor 6xy is omitted (as it was done in Eq. (j26|l ) since all 
the amplitudes are diagonal in the proton's spin. 

II.5 THE BFKL LADDER 

As we shall see below, the characteristic transverse momenta in our gluon loop are 
large, so the estimate is sufficient for our purposes. Substituting the nucleon 
impact factor into Eq. ()12|) one obtains: 

47r2 kj_{r + kf 

The final integration over k^ reveals the logarithmic dependence of the photon impact 
factor lnr2/l]2 _^ hi\t\/Q'^: 

cPki_ f 1 1 \ 

kl{k^ + r^)2 \{k^ + ari_f + fi2 ^2^2 ^ ^2 j 

(21" T^ + 'l) 



V-'i^'iEe;Writ)l'^J^. (42) 



(i2fc_L / (fc_L + ar_L)'^ (-r_L + ar_L)^\ 



A;i(A;x + r±)2 V(A;x + arx)2 + fi2 ^2^2 ^ ^]2 



TT rl^ l]2_^2^2 ^2 ^ ^2^2 ^2 _ ^2^2 f]2^^Vi 

fi2-aar2 ;f Ifi2 + ct2r2 ^ WT^h^ ^ fi2 

fi2 - a^rl rl o?r\ o?r\ > 

^?]2^"^ 1]2 + " ?]^T^ ~ l^^T^^''-' ' ^ ^ 

/■ d2fcx / {k^ + ar^)2 ^2^2 \ 

7 A;2(A;^ + r^)2l (fc^ + ar^)2 + (]2 a2r2 + fi2 ) ^O' ^^^^ 
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(fk± ( (k± + ar±)'^{k± + ar±Y a^rj_{-r± + ar±) 



M 2 

irr, r 



-fi^*/i + ^(ln^ + ry), (46) 



where 



T] = ^ + -fE + Inn + \nn'^, (47) 
R^^Ru ^ \9^uR\ (48) 

and is defined in Eq. 

For the case of transverse photon polarizations one gets 

V^^ = §f^)'(Ee?)^rW (49) 
q 



X \ TT 



^ In^ |[ + 2) - (e, 6')x + ^(e, r)^ie', r)^ + 0(^ 



and when the incoming photon has longitudinal polarization we have 



'iW^_5,„^i; + li-^ + 0(-l)). (50) 



The longitudinal amplitude (j^Uj) is twist-suppressed as y\t\/Q'^ in comparison to the 
transverse amplitude (j5(jp (as it should, due to the fact that t ^ corresponds to a 
real incoming photon). 

Since the integral over k± fl42|) converges at k±_ ~ Q, the region k± ~ M, where 
we do not know the nucleon impact factor, contributes to the terms ~ (9(|t|/(5^) 
which we neglect. 

In the next-to-leading order (NLO) in perturbation theory the most important 
diagrams are those of the type shown in Fig. El Actually, this diagram gives the total 
contribution in LLA if one replaces the three-gluon vertex in Fig. El by the effective 
Lipatov's vertex [22]. 

Calculation of this diagrams in the leading log approximation yields: 

£k^£k'^ I{k^,r^) lNik'^,rj_) 
47r2 47r2 A;i(r + A;)i^^'^'^\fcl)2(r + fc')i' 
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FIG. 9: Typical diagram in the next-to-leading order in perturbation theory. 



where K{k±, k'j_,r±) is the BFKL kernel [30] 

^k'-{k' + r) k'-{k'-k) {k' - k) ■ {k' + r) 



K{k',k,r) = 2 
4 



k'^{k' + ry k'^{k'-ky {k' -ky{k' + ry 

k' ■ {k' - k) {k' - k) ■ {k' + r) 



+ 



(52) 



{k' - ky[k'^ + {k' - ky] {k' ~ ky[{k' - ky + {k' + ry]_ 

As we shall see below, the integral over k'j_ converges at \k'j_\ ^ M so we can again 
use the approximation (j^Bj) for the nucleon impact factor. One obtains 



J d'k'^K{k^,k'^,r^) — 



7rFf+"(t) (ln::^ + ln 



k] 



{k'^y{r + k')l 
and therefore the amplitude (j5T| takes the form 



{k + rf_ 



{k + r)[ 



Finally, the integration over k yields: 



V 



2 fa 



X \ TV 



5 - iT.4)Fr{t) 



3a 



lu- 
ll X 



((.,.')x(iln'g + 21„g-2 + C(3)) + (^(., 



r±(e,rh - (e, e 



(53) 



(54) 



(55) 
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where the accuracy is 0(1/ Inx) and ({3) is a Riemann zeta function. 

In the next order in the BFKL approximation (see Fig. ITUj) it is still possible to 
obtain the DVCS amplitude in the explicit form. It is possible to write down the 




FIG. 10: Typical diagram in the next-to-next-to-leading order in perturbation the- 
ory. 

result of the summation of the BFKL ladder in the form of a Mellin integral over 
complex momenta using the Lipatov's conformal eigenfunctions of the BFKL equa- 
tion in the coordinate space. Unfortunately, we were not able to perform explicitly 
the integration of the Lipatov's eigenfunctions with impact factors and without it the 
Mellin representation of the DVCS amplitude is useless for practical applications. 
The amplitude in the next-to-next-to-leading order is: 

^ = -(Eej(6a.ln-j / ^^^^ ^ (56) 

Once again, if we use the fact that the integral over k'^ converges at \k'^\ ^ M we 
can approximate the nucleon impact factor by Eq. (^Dj), and obtain: 

r<Pk'^ fd'k'i 1 Mk'lr^) 

J 1^7 -^i^(fcx,A.^,r^) ,i^(fc^,A;^,r^) 



(A;l)2(r + k')l" ■ {k'iy{r + k")l 

ik\r {k' + r)l \ 
{k')l{r + k')l rl rl J 



1 M'fcl K{k^,k'^,r^) ( 



22 



(k + r) 



The resulting integration over k± yields: 



(57) 



V 



9 

X \ TT 

2(C(3) 



X 



■V24 



In^ 



+ In' 



2 In 



1) + 1.46) + ( -^(e,r)^(e',r) 



(5^ 



As it was mentioned, we are unable yet to obtain the explicit expressions for the 
amplitude in higher orders in perturbation theory. It turns out, however, that for 
HERA energies the achieved accuracy is reasonably good. The estimation of the next 
term gives ~ 30% of the answer at not too low x (see the discussion in next section). 

In the leading logarithmic approximation it is impossible to distinguish between 
«s(Q) and a;s(y^) - to this end one needs to use the NLO BFKL approximation [31] 
(see also [32]) which is beyond the scope of this paper. 

The final result for the DVCS amplitude with transversely polarized photons is: 



V 



(59) 



2 fasiQY 



X 



TT 



where 



vix,Qyt) 



- In^ ^ + 2 + 



{e,e')^v+^{e,r)^{e',r)^-{e,e')Av' 



t 



2 V TT X 



24 



In^ 



Q 



+ In^ ^ + 2(C(3)-1) In 



1.46 



v'ix, QVt) = 1 + In 1 + i f In 

X 2 \ TT X , 



TT 



(60) 
(61) 



Note that the spin-dependent part ~ v' does not contain any ln(5^/|t| and is, 
hence, much smaller than the spin-independent part ~ v. For the longitudinal polar- 



ization the amplitude is twist-suppressed as ~ y so we have not calculated 

any terms beyond Eq. ()50|) . In the numerical analysis carried out in the next sections 
only the spin- independent part of the amplitude is kept: 

^\Y.^^e'^y = I (^)^Ee?)^^"w^(^.Q^^)• (62) 
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The expressions above give us the imaginary part of the DVCS amphtude. For 
the calculation of the DVCS cross section one needs to know also the real part of 
this amplitude, which can be estimated via the dispersion relation. For the positive- 
signature amplitude H± (= \j2^±^'j_H) we get [33] (see also [19]) 



which amounts to the substitution 

In s — i> 



^tan(.^) '^H^is), 



ln(— s — ie) + Ins) 



in our amplitude (jU^ . Thus, the real part is: 



R = -^Hj 

n 



2 

X \ n 



{j:el)Fnt)r{x,Q\t), 



(63) 



(64) 



(65) 



TT 



r{x,Q\t) = - 



3a. 



TT 



_ln^:^ + 21n:^-2 + C(3) 



3a. 



TT 



In 



X V24 



In4^ + In2^ + 2(C(3) 

\t\ \t\ 



l)ln^ + 1.46 



(66) 



II.6 DVCS CROSS SECTION 



It is instructive to compare the DVCS amplitude V given by Eq. © with the corre- 
sponding amplitude for the forward 7* scattering 



H 



te^e' / dze-''^-^{p\T{f{z)fm\p). 



(67) 



The imaginary part of this amplitude is the total cross section for deep inelastic 
scattering (DIS) 

1 



-^H = W 



TT 



(68) 



For example, W averaged over the transverse polarizations of the photons is: 



(69) 



(at the leading twist level we have the Callan-Gross relation F2 = 2xFi). We will 
compare the imaginary part of the DVCS amplitude V± given by Eq. (jU^ to the result 



24 



for W± calculated with the same accuracy. (The notation W±{x) is used rather than 
Fi{x) in order to avoid confusion with Fi{t)). 

Similarly to the DVCS case, the DIS amplitude has the form (cf. Eqs. (jl^ . (fKT|l . 
and dMl)): 

2g^s ON rd^ki 1 



1 ^ 1 



, K{k^,k'^,0) f,2,nK{k'^,k'i,0) 



(70) 



where l±{k±,0) is the virtual photon impact factor averaged over the transverse 
polarizations [34]: 



/±(fc±,o) 



1 /•! da /■! da' kl{l - 2aa){l - 2a' a') 



(71) 



2 7o 27r 7o 27r a'a'(A;2 + 

The nucleon impact factor /7v(A;^,0) cannot be calculated in perturbation theory 
since it is determined by the large-scale nucleon dynamics. However, we know the 
asymptotics at large k^^ M 



lN{k±,0) 



fc2 »A/2 



Ff+"(0) = 1. 



(72) 



Also, at lN{k±, 0) — * at A; ^ due to the gauge invariance. It seems reasonable to 
model this impact factor by the simple formula: 

k', 



iNik^.O) 



k^, +M^' 



(73) 



which has the correct behavior both at large and small k±. With this model, the 
DIS amplitude (fTUj) takes the form: 

2 



w± = — = — 
2x 3x 



asiQ) 



3a, 



IT 



TT 
")2 



Q , 7, Q 



77 



7 

M2 ' 12 



- In^ + - In H 

2 M2 6 M2 

77 , 



(74) 



M2 ' 18 



27 



2 V vr X 



24 M2 36 M2 36 



131 ... ^\, 1396 TT^ 14,, ' 
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Note that the coefficients in front of leading logs of Q^, determined by the anomalous 
dimensions of twist-2 operators, coincide up to the factor 2/3. The graph of the 
model ()74|1 versus the experimental data is presented in Fig. ^2 for = lOGeV^ 
and g2 = 35GeV2 (we take Ee^ = f ). 



F2 



Q**2=10 



F2 




Q**2=35 




3 -2 -1 



Log X 



FIG. 11: F2{x) from Eq. (fflj) versus experimental data at Q"^ = lOGeV^ and Q"^ 
35GeV2. 



For DIS it is possible to write down the total BFKL sum as a Mellin integral and 
unlike DVCS, the integrals of impact factors with the BFKL eigenf unctions (fc^)"^"^*'^ 
can be calculated explicitly. Eqs. (fflj) and (f73j) correspond to the expansion of this 
explicit expression in powers of aslnx. 

For example, the next term in BFKL series ()74|) has the form: 



4 

3x 



1 



TT ) 
,5 



(E e: 

flavors 

7 



■In" ^ + 
120 M2 144 



In' 



, 1 
m — 



(75) 



M2 



77 , 3 

+ — In^ ^ + 



+ ( 



,1396 



^ + ^C(3)).n|; 



108 M2 
4736 Itx^ 



243 



g„+fc(3) + 6C(5) 



The ratio of this (a^ Inx)'^ term to the sum of the first three ones ()74j) is presented 
in Fig. ^Jfor = lOGeV^ and = 35GeV2. From these graphs we see that the 
sum of the first three terms gives the reliable estimate of the DIS amplitude at not 
too low X. It is expected that the same will also be true for the DVCS amplitude. 

At very small x ~ 10^'^ — 10~^ the full BFKL result for F2 in our model is 
growing more rapidly than Fig. El On the other hand if one takes into account 
the NLO BFKL corrections [31,32] the result for F2 at very small x goes well under 
the experimental points. This indicates, that at such x we need to unitarize the 
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BFKL pomeron, which is currently an unsolved problem. (The best hope is to find 
the effective action for the BFKL pomeron (see e.g. [35,36])). On the contrary, at 
"intermediate" x ~ 0.1 — 0.001, we see from Fig. ^Jthat, since the corrections almost 
cancel each other, it makes sense to take into account only a few first terms in BFKL 
series. 

It is instructive to compare the t-dependence of our DVCS amplitude (jUnj) with 
the model used in the paper [19] 

V^{x,t,Q') = }-F,{x,Q')e'"\ (76) 
V2(x,t,Q2) = ^Fi(x,g^) (77) 

where R ^ 0.5 for our energies. (Literally, the model used in ref. [19] corresponds to 
Vi but it is more natural to approximate the t - dependence by the dipole formula 
[37]). 

The comparison is shown in Fig. Elfor = lOGeV^, = 35GeV^ and a;=0.01, 
a;=0.001. 

In order to estimate the cross section for DVCS at HERA kinematics (Q^ > 6GeV^ 
and X < 10"^) we will use formulae from Ref. [19] (see also Ref. [38]) with the trivial 
substitution -^F2{x)R~^e^^^'^ Vj_(a;, Q^, t). The expressions for the DVCS cross 
section, the quantum electrodynamics (QED) Compton (Bethe-Heitler) cross section, 
and the interference term have the form {y = 1 — y): 



dxdydtd(f)r Q^y 
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dxdydtdcpr 



(1 + w^) 

T2a^ ^ t)i^f (t) cos 0,. 



(79) 
(80) 



dxdydtdcpr Q^\/y\t 
The expression for the interference term from ref. [19] is corrected by factor 2 [37] 



x=0.01,Q**2=10 



x=0.001,Q**2=10 
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FIG. 13: The ratio Vi/V± (lower curve) and V2/V± (upper curve). 



[39]. Here y 



{E and E' are the incident and scattered electron energies, 



respectively, as defined in the proton rest frame) and (j)r = (pe + <Pn where 0Ar is the 
azimuthal angle between the plane defined by 7* and the final state proton and the 
X — z plane and 0e is the azimuthal angle between the plane defined by the initial 
and final state electron and x ~ z plane (see Ref. [19]). As mentioned above, we 
approximate the Dirac and Pauli form factors of the proton by the dipole formulas 

At first let us discuss the relative weight of the above cross sections. We start 
with the asymmetry defined in ref. [40] 



A 



/2 



Cd<Prda^^^ 



^1) 
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Asymmetry at x=0 . 01 , Q**2=10 Asymmetry at x=0 . 001 , Q**2=l 




FIG. 14: Asymmetry versus y = 0.1 — 0.6 and |t| = 1 — 5 GeV^. 



where 

^ dcT^^^s + da^^D^ + da'""^. (82) 

The asymmetry shows the relative importance of the interference term, which is 
proportional to the real part of the DVCS amplitude. In our approximation the 
asymmetry is 

%/i(Ee^) (f)'(l + 2.8J^)r 

A{y,t) = ^—^ , . (83) 

47r2(E einv^ + r^) [^) [l + ^) + {} + 7.84 

The plots of asymmetry versus y and \t\ are given by Fig. 

Second, we define the ratio of the DVCS and Bethe-Heitler cross sections [19] 

= 'd^ = „2n2 /. . ^ ^ . 1,1 X • (84) 

This ratio is presented on Fig. ^1 









^ + 4M2 ) 


1^(1 + 7.. 
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We see that there is a sharp dependence on y; at y > 0.2 the DVCS part is 
neghgible in comparison to Bethe-Heitler background whereas at y < 0.05 the QEDC 
background is small in comparison to DVCS. 

Finally let us estimate the relative weight of the DVCS signal (starting from 
\t\ = 1 GeV^) as compared to the DIS background. We define (cf. ref. [19]) 

^ ^ a{Y + P ^ 1 + P) _ 
0"(7* + p — > 7* + p) 



Q2F2(x,Q2) 



(f) (Eeg lfdt{Fr{t)) {v\x,Qyt)+r\x,Qyt).) (85) 



-5 



At g2 = lOCeV^ we find = 1.56 x 10"^ for x = 0.01 and Ry = 2.36 x 10 
for X = 0.001, while for = 35GeV^ we find R^ = 0.62 x 10"^ for x = 0.01 and 
R^ = 0.71 X 10"^ for X = 0.001. 

The expressions (l7H jl - (|Hn|) are correct if ^ |t| up to 0{-^) accuracy with the 

notable exception of the correction 0(y'^) coming from the expansion of electron 
propagator in the u-channel of the Bethe-Heitler amphtude. As suggested in ref. [39], 
at intermediate t one can keep the propagator in unexpanded form (and expand the 
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FIG. 16: Asymmetry with the correction factor (j86|) . 



rest of the amphtude, as we have done above). This amounts to the replacement 



y^y 



Q^y 



\t\y ^ 



,(2-2/) [\t\ 



cos- 



1^1 9 

+ cos' 

Q 



6) 



Vy 

in the numerator in Eqs. (f7^ and (IHUj) (see ref. [38]). 

The resulting asymmetry ()81|1 is presented in Fig. El We see that the correction 
factor (j86|) crucially changes the behavior of the asymmetry due to the fact that it 
restores the azimuthal dependence of the QEDC amplitude which was not taken into 
account in Eqs. ()78II80|) . In order to find asymmetry at these and t with greater 
accuracy one should take into account other twist-4 contributions as well. On the 
contrary, the ratio D{x,Q'^/t) does not change much (see Fig. ITTj) so we hope that 
our leading-twist results for the ratio presented in Fig. [Tulare reliable. 



II.7 CONCLUSION 



The DVCS in the kinematical region (0)) is probably the best place to test the mo- 
mentum transfer dependence of the BFKL pomeron. Without this dependence, the 
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model (|77|) would be exact, hence the upper curves in Fig. indicate how impor- 
tant is the t-dynamics of the pomeron. We see that the t-dependence of the BFKL 
pomeron changes the cross section at t > 2GeV^ by orders of magnitude and therefore 
it should be possible to detect it [41]. 

The pQCD calculation of the DVCS amplitude in the region (@)) is in a sense 
even more reliable than the calculation of usual DIS amplitudes since it does not 
rely on the models for nucleon parton distribution. Indeed, all the non-perturbative 
nucleon input is contained in the Dirac form factor of the nucleon, which is known to 
a reasonable accuracy. There are, of course, the non-perturbative corrections to the 
BFKL pomeron itself. It is not clear how to take them into account at this moment. 
Of course, any reasonable models of nucleon parton distributions such as ()32|) should 
reproduce the form factor after integration over X. 

Finally, let us discuss uncertainties in our approximation and possible ways to 
improve it. One obvious improvement would be to calculate (at least numerically) 
the next ~ (oslnx)^ term in the BFKL series for the DVCS amplitude. Hopefully, 
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it will be as small as the corresponding calculation of the DIS amplitude suggests. 
Besides, there are non-perturbative corrections to the BFKL pomeron which are 
mentioned above. These non-perturbative corrections correspond to the situation 
like the "aligned jet model" when one of the two gluons in Fig. El is soft and all 
the momentum transfers through the other gluon. It is not clear how to take these 
corrections into account, but one should expect them to be smaller than the corre- 
sponding corrections to F2{x) which come from two non-perturbative gluons in Fig. 
El (in other words, from the "soft pomeron" contribution to F2{x)). 

The biggest uncertainty in our calculation is the argument of coupling constant 
Qs which we take to be Q^. As mentioned above, it is not possible to fix the argument 
of as in the LLA, so one could have used Q;s(|t|) instead. We hope to overcome this 
difficulty by using the results of NLO BFKL in our future work. 
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CHAPTER III 
GENERALIZED PARTON DISTRIBUTIONS 

III.l INTRODUCTION 

Let me start with the historical overview of the problem. The first information 
about quark distributions inside the nucleon was received from the experiments on 
deep inelastic scattering: 

eik) + N{p)^e{k')+X{pn), (87) 

where X{pn) are possible final hadronic states with the 4- momentum p„. 
The scattering amplitude for the n-channel is given by: 

T„ = eMk',X'h^u{k,X)^{n \ J,^^(0) | p,a), (88) 

where e is electron's charge, k{k') is the initial(final) momentum of the electron, A(A') 
is electron's initial(final) polarization, q is the momentum of the virtual photon, that 
hits the nucleon (g^ = {k — k'Y = ^4£'£''sin^| < 0, = —q^), 9 is electron's 
scattering angle and is hadronic electromagnetic current. 

The differential cross section for this scattering for unpolarized particles has the 
form: 



111 d^k 



3 IJ n 



I ^; I 2M 2^(271)32 



d^Pi 



(27r)32p,o 



^ |7;n27r)^(5^(p + A;-A;'-p„), 

(89) 

where Pn = EILiPi- 

If we sum over all possible final hadronic states, we will get the inclusive cross 
section 

where a — e^/An - fine structure constant. Lepton tensor Ifj,^, has the following form: 

l^, = ^Tr[k%h.] = 2 (kX + k'X + y^/^-) ' (91) 
and the hadronic one W^^ has the form: 

= 4^ E / ^ I {J'e.m{^)J:.mm I P, ^) ■ (92) 
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Prom the current conservation dfj,J^^ — one gets the following equations: 

q^^W^'' = QuW^'' = 0, (93) 

which brings us to the requirements that hadronic tensor W^''' should satisfy: it 
should be a Lorentz invariant tensor of rank two and it should depend on momenta 
Pi^ and (the only momenta we have in this case) , which brings us to the following 
representation of the hadronic tensor: 



-Wi (g^, - ^ ) + Ip^ - (p. - 



(94) 



where Wi^2 are the first Lorentz invariant structure functions of the nucleon target 
that had ever been introduced. Those structure functions depend on two variables 
q'^ andu^E- E' . 

In the case, when the final hadronic state X{pn) is also a nucleon (the case of 
the elastic electron- nucleon scattering), the matrix element for the electromagnetic 
current has the form: 

1 



{Nip') I Je^„(0) I Nip)) = nip') 



^^Fr{q')+ta,,q''—F2{q' 



u{p), (95) 



with q — p — p' . For the proton Ff (0) = 1 is equal to the electric charge and 
F^O) = 1.79 is equal to the anomalous magnetic moment. For the neutron those 
number are -F"(0) — and -F2"(0) — —1.91, respectively. 

Therefore, the differential cross section measurements for the elastic electron- 
nucleon scattering give us information about electrical and magnetic form factors. 
If the nucleon had no internal structure, those form factors were independent on q^. 
The experimental fact that they have momentum transfer dependence opened a new 
era in nucleon structure investigations. 

The new parton model was created by Feynman, which described a nucleon as a 
particle consisting of three partons (quarks) interacting by gluon exchange. Lepton- 
nucleon scattering has been transformed to the scattering of virtual photon off one of 
the partons. Accordingly, the structure functions Wi^2i'Pil) describing the response 
of the nucleon with momentum p to the momentum transfer q were replaced by 
structure functions Fi^2{x): 

MWi (p, q) ^ Fi (x) = ^n{x) , (96) 
uW2{p,q)^F2{x)^xn{x), (97) 
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FIG. 18: Forward quark-parton density (left) and helicity (right) distributions used 
in our calculations for u (solid) and d (dashed) quark. 



where n{x) represents the probability to find the parton (quark) inside the nucleon 
carrying the fraction x of initial longitudinal nucleon momenta xp within the range 
from X to X + dx. The corresponding function for the antiquark is denoted as n{x). 

Thus, through the experiments on deep inelastic lepton-hadron scattering one can 
extract information on forward quark distribution functions (Fig. ITHj) . 

In the following, we introduce the kinematical variables entering the DVCS pro- 
cess. We then discuss the DVCS amplitude in the leading power in Q (twist-2 
accuracy), and show its dependence upon the general parton distribution (GPD) 
functions. 

III.2 DVCS IN TWIST-2 APPROXIMATION 

As was already mentioned in the first chapter, the deeply virtual Compton scattering 
process 

^*(g) + iV(p)^7(g') + iV(p'), (98) 

in the limit of vanishing momentum transfer t = {p—p'Y ^"^^ large virtuality ^ oo 
offers us much more complex insight on parton dynamics in the nucleon. In the lowest 
approximation (twist-2) the structure functions F{x,y), called double distributions 
[1,2], carry information on the quark distributions with different fraction x of the 
initial nucleon momentum and different fraction y of the momentum transfer r = 
p' — p carried by the quark. This nucleon structure information can be parametrized 
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FIG. 19: Double distributions F{x,y;t = 0) for u (left) and d (right) quark. 

in terms of four structure functions. Those functions look like usual distribution 
functions with respect to x and like distribution amplitudes with respect to y. 

III. 2.1 Kinematical variables and distribution functions 

Since q'^ = (g' - is the final real photon four-momenta), it is natural to use it as 
one of the basic Sudakov light-cone 4-vectors. Another basic light-cone vector we 
use is p, which is also understandable, since = can be neglected compared to 
the virtuality = — of the initial photon and the energy invariant p ■ q = Mv. 
In this limit (p^ = and = 0), the simple equation p'^ = (p + — requires 
p ■ r = which can be satisfied in the light-cone basis only if they are proportional 
to each other r = Cp. 

Despite the proportionality between p and r, they specify the momentum flow in 
two different channels, s and t respectively. The coefficient of their proportionality C 
is the parameter characterizing the "asymmetry" or the "skewedness" of the matrix 
elements. 

The total momentum fraction of the quark can either be positive or negative. 
Since positive momentum fraction corresponds to quark and the negative one - to 
antiquark, there are three major regions that could be found. Two of them, namely 
when we have two active quarks or antiquarks, are the regions accessible in inclu- 
sive processes like DIS and giving us the usual forward parton distributions. The 
third region represents the one-quark-one-antiquark case, which cannot be measured 
through DIS, since it vanishes in the forward limit A ^ 0. In the limit r = we 
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come up with the following reduction formulas for the double distributions F{x,y): 

rl—x rl 

/ F{x,y;t = Q)\^>ody = n{x); / F{x,y;t = Q)\^<ody = -n{-x). (99) 

Jo J-x 



/ F{x,y;t^O)\^>ody^n{x); / F{x,y;t ^ 0)\^<ody ^ n{-x). (100) 

Jo J-x 

Positive-x and negative-x components of the double distributions can be treated 
as non-forward generalizations of quark and antiquark densities, respectively: 
F{x,y;t)\^>o {-F{-x, 1 -y;t)\^<o)- 

The reduction formulas and interpretation of the x-variable suggests, that the 
profile of F{x,y) in x-direction is basically determined by the shape of n{x). The 




FIG. 20: "Handbag" diagrams for DVCS. 

profile in the y-direction at the same time characterizes the spread of the parton 
momentum induced by the momentum transfer r. This brings us to the following 
relation: 

F{x, y) = h{x, y)n{x), (101) 
where h{x,y) is an even function of y and 

"l-|rr| 
-l+\x\ 

Integrating each particular double distribution over y gives the non-forward par- 
ton distributions: 

nl rl-x 



[' h{x,y)dy^l. (102) 

J-l+\x\ 



Tc{X)=[ dxj ^ S{x + Cy-X)F{x,y)dy^ 
Jo Jo 

e{X > C) / Fix - Cy, y)dy + e{X < C) / F{X - Cy, y)dy, (103) 
Jo Jo 
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where X = x + (y is the total momentum fraction and the notation = 1 — is 
used. 



III.2.2 Twist-2 DVCS amplitude 

The leading [1,3,42] handbag contribution fFig. EUj) to the DVCS amplitude can be 
represented as [14] 



2{pq 

X 



1 



+ 



ief'^^f^p^q'^ 

{pq') 

u{p')q'lMp)Tm + ^-^u{p')lMp)Tm 



(104) 



X 



The invariant amplitudes T"((^) can be calculated using model [43] for the non- 
forward quark parton distributions with t = 0: 

»i 



1 



(105) 



_X-C + ie X-ie_ 

Note, that because the non-forward distributions are real, the imaginary part of 
T"(C) comes only from the singularities of the expression in the square brackets. 
Since all non-forward distributions vanish at X = 0, only the first term in the square 
brackets generates the imaginary part: 

The real part of the invariant amplitudes T"-{C) (Fig. I^T|) is given by the "principle 
value integration" in Eq. ()105|l 

dX 



The non-forward distributions JF"(X) are expressed through the double distributions 
by Eq.dMI). 

In our calculations we use the model [43] for the double distribution given by the 
factorized ansatz 



^a{x,y) = fa{x) ^^^/^\ h{x) = f h{x,y)dy, 
h[x) Jo 



(108) 
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where fa{x) is the forward quark - parton distribution and h{x,y) is the profile 
function. The reahstic profile for the double quark distributions, both valence and 
sea, is close to the "asymptotic" form 

has{x,y) = 6yil-x-y). (109) 

The resulting curves for the real and imaginary parts of the invariant amplitude 
T{0 = EaelT^iC) are shown in Fig. 121 

III.2.3 Twist-2 DVCS cross section 

The DVCS amplitude, as was mentioned before, may be observed in the exclusive 
lepton-nucleon scattering. In the 7*(g) + N{p) l{<i') + N{p') reaction, the final 
photon can be emitted either by the proton or by the lepton. The three relevant 
diagrams are shown in Figl^ The blob with two photon legs in diagram FigQHh 
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FIG. 22: Real (solid) and imaginary (dashed) parts of the DVCS invariant ampli- 
tudes: Tp (left top), Tk (right top), Tq (left bottom), and Tp (right bottom). 

stands for the DVCS nucleon amplitude (i.e. the amplitude T^vcs (in Eq. (jllOp ) 
of scattering of the virtual photon on the nucleon, with the real photon in the final 
state, whereas the one-photon blob in diagrams Fig]23L'. c represents the nucleon 
electromagnetic form factor. 

We will refer to diagram FigQHh as the DVCS part of the amplitude of the 
electron-nucleon scattering. Figsl23b. c give together the Bethe-Heitler part. 

The invariant cross section of the '~^*{q)+N{p) l{(l')+N{p') reaction differential 
with respect to Q^, xb, t, and out-of- plane angle (0 = 0° corresponds to the 
situation where the real photon is emitted in the same half plane as the leptons) is 
given by : 

da _ 1 xbV^ f AM^xiy 

dQ^dxBdtd<^ ~ 32 \ ^ J 



Tbh + Ti 



DVCS 



(110) 
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FIG. 23: a) the DVCS part of the amphtude; b),c) the Bethe-Heitler part. 



where M is the nucleon mass, y = {p ■ q)/{p ■ k), and k is the initial lepton four- 
momentum. 

The process where the photon is emitted from the initial or final lepton is referred 
to as the Bethe-Heitler (BH) process (amplitude Tbh in Eq- piO|) ^- and can be 
calculated exactly: 

,6 



T (111) 



(4l/V(l + 8/5' + t' + 4/5(1 + r))((2Fi(i) + F^it)^ + 2(2/i/r - l)F^{t)) 



+r(l + y'- XBy{l + 2/5(1 + y) + yT)){A^Ft{t) - tF. 



where we introduced the dimensionless variables (3 = {k ■ q')/Q'^, t = t/Q^, and 
fi = M'^/Q'^. The 0-dependence is hidden in /5 and it alone determines the 0- 
dependence of the Bethe-Heitler amplitude Tbh- 

Light particles such as electrons radiate much more than the heavy proton. There- 
fore the BH process generally dominates the DVCS amplitude, especially at small 
t. The best way to extract the information on general parton distribution functions 
and still be in an accessible region for the experimentalists is to measure the in- 
terference term. We could benefit from the situation if we consider the difference in 
cross section for electrons with opposite helicities. This way, in the difference of cross 
sections a^^ — cTg^, the BH (whose amplitude is purely real) drops out. The DVCS 
contribution is strongly suppressed and the term dominating the cross section would 
be the BH-DVCS interference [44] term 



rpBHrpDVCS* 



(112) 
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where the Bethe-Heitler process will project out the imaginary part of DVCS ampli- 
tude magnifying it with its own full magnitude. 

The imaginary part of the DVCS amplitudes is directly proportional to the GPDs, 
evaluated along the line x = ^, and measures in this way the 'envelope functions' 
H{^,^,t), E{^,^,t), Hi^,^,t), and (Eq. »)• 

Figures El through 1221 show the results of our calculations for the energies and 
kinematics accessible at Jefferson Lab. The values of transferred momentum squared 
run from t = —.1 to t = —.5 limited by 

t^,„ = -4M2eV(l-e') (113) 

for each of kinematical sets. 

One can see, that the interference term is larger at lower momentum transfer 
t even though it is responsible for a smaller percentage of the total cross section. 
As the initial energy of the electron E^. goes up, the interference cross section value 
drops down by one order. The ratio of the interference term to the total cross section, 
nevertheless, stays the same: about 40 to 60 percents. 

III.3 TWIST-3 CORRECTIONS FOR THE DVCS PROCESS 

The twist-2 DVCS amplitude ()104|) is independent of Q. This means it has a scaling 
behavior. However, this amplitude is not complete beyond the leading order in 
Q. Although it is exactly gauge invariant q^T^^^ = with respect to the virtual 
photon, the electromagnetic gauge invariance is violated by the real photon except 
in the forward direction t = 0. This violation of gauge invariance is a higher twist 
(twist-3) effect compared to the leading order term T'^'^. Since the product of final 
photon 4-momentum and the DVCS amplitude at twist-2 accuracy is proportional 
to the transverse component of the transferred momentum A^, an improved DVCS 
amplitude has been used [45,46] to restore the gauge invariance in the nonforward 
direction: 

pu 

^DVCS = ^DVCS-2 + JJdT^ (^-l)a ^DVCS-2 • (H^) 

It creates a correction term of higher order in Q to the twist-2 DVCS amplitude 

-'dVCS-2- 

It is necessary to have an estimate of the effect of power suppressed (higher twist) 
contributions to those observables in order to be able to extract the twist-2 GPDs 
from DVCS observables at accessible values of the hard scale Q. 




FIG. 24: Cross section calculated at twist-2 level for Eg — 4.25GeV, xb — -2, and 
= 1.3GeV^ for the different values of momentum transfer: t — —.1 - solid curve, 
t — —.2 - dotted curve, t — —.3 - dashed curve, and t — —.5 - dash-dotted curve. 
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FIG. 29: The same for Ee = 6GeV, xb = -3, and = l.SGeV^. 
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FIG. 30: The same for = 6GeV, xb = -3, and ^ 2GeV2. 
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FIG. 31: The same for = 6GeV, xb = -3, and ^ 2.5GeV2. 
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The first power correction to the DVCS amphtude of order 0{1/Q) is called 
twist-3. 

The twist-3 corrections to the DVCS amplitude, which have been derived and 
calculated recently by several groups [47-56] using different approaches. In our study 
we used the notations and definitions from [57]. The following section is a short 
overview of these definitions. 

III. 3.1 Kinematical variables of the DVCS process 

Let us start with the comparison of the kinematical variables used for the description 
of twist-3 amplitude with the set of variables we used in the previous section. There 
is only a slight difference, but it is worth mentioning. 

There are four generalized structure functions used for the description of DVCS 
amplitude. They are denoted by H, E, H and E and depend on three variables. The 
light-cone momentum fraction x is defined by = xP~^, where k is the quark loop 
momentum and P is the average nucleon momentum (P = {p + p')/2, where p {p') 
are the initial (final) nucleon four-momenta respectively). In the previous chapter 
we have used the initial momentum of the hadron p as one of the basic light-cone 
momenta and up to twist-3 accuracy those momentum coincide as will be shown 
later. The skewedness variable ^ is defined by A+ = —2^P^, where A = p' — p is 
the same as r in the previous chapter, the overall momentum transfer in the process, 
and where 2^ — xb/2) in the Bjorken limit. The third variable entering the 

GPDs is given by the Mandelstam invariant t = A^, exactly as it was before. 

Using momenta P and A instead ofp = P — A/2 and p' = P + A/2 we can 
define the active quark momentum as k^ — A+/2 = (x + ^)P^ before the virtual 
photon impact and k~^ + A+/2 = (x — ^)P^ afterwards. The light-cone basis is 
chosen along the positive and negative z-direction (both and P^ are coUinear 
along the z-axis and have opposite direction): p^ = P+/-\/2(l, 0, 0, 1) and n'^ == 
1/P+ ■ l/\/2(l, 0, 0, — 1), satisfying p-p = 0, n-n = 0, p-n = l and using the 
notation P+ = P ■ n. 

In this frame, the physical momenta entering the DVCS process fj98|) have the 
following decomposition [45] : 



1 777^ 



(115) 




(116) 



53 



± 5 



(117) 
(118) 



where is the perpendicular component of the momentum transfer A ( i.e. p-A_ 
n ■ A_|_ = ), and where the variables fh?, ^' and ^ are given by 



A^ 

^2 ^ - — 
4 



2^' = -q' 



P-q 



1 + 



2^ = -A+ = 2^ 



g^-A 



Bj Xb 



Bj Xb 



1 XB ' 
2 



To twist-3 accuracy, Eqs. ()115II118|) reduce to 



1 — £B 
2 



P 
A 



(119) 
(120) 
(121) 



(122) 



-2^P + A^, 
q = -2^P+^n, 

q = -^n-A^. 

Thus, the final (real) photon 4-momentum squared g'^— > A\ suggests the twist- 
3 accuracy t = A^ ^ A^ 0. 



III.3.2 Twist-3 DVCS amplitude 

As a starting point we use the DVCS amplitude on the nucleon to the order 0{1/Q) 
derived explicitly in a parton model approach [48] and in a light-cone expansion 
framework [49]: 



H 



1 



2 7-1 



(P-q) 



dx 



-9 



{P ■ q') 



- {P-q' 

pu^k 



n%(x,OC+(x,0 (123) 



(P ■ q') 



{^fc(x,0 C+(x,0 
-i{eu,)^P{x,0 C-{x,i)] 
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where the functions and Tfj, are given by: 

^^(^,0 = ^^N{p')N{p)E{x,i) - ^N{p')hN{p){H + E){x,0 

+ £^duG^{u,0W4x,u,0 + te^^k j'^duG\u,0W4x,u,0, (124) 

^,.(^,0 = ^N{p')^,N{p)E{x,0 - ^N{p')h^,N{p)H{x,i) 

+ f_^duG^{u,OW+{x,u,0 + ie^^.k j'^du G''{u,^)W_{x,u,0 ■ (125) 
The following notations are used : 

^•(u, = N(p')jiN(p)(H + E)(u, 



. du 



^N{p')hN{p) 



. du d^. 



{H + E){u,0, (126) 



G^{u,^)^N{p')^^^,N{p)H{u,^) 



A'' 



d d 



^N(p')h^,N{p) 



u— + 

du 



H{u,0. (127) 



The functions W±{x,u,^) arc called [58] Wandzura-Wilczek kernels. They were in- 
troduced in Ref. [49,51] and are defined as : 

/ 1 f/i/ ^.(^iu>x) ^, ^.9(u<x)'] 

W±{x,u,0 = ::\0{x>O- -^-e{x<0- ^ 

21 u — t u — t ) 



±-< dix > -£)— - 9(x < -£)— >. 

2l ^ u + ^ ^ u + ^ J 



(128) 



We also introduce the metric and totally antisymmetric tensors in the two dimen- 
sional transverse plane (£0123 = +1) : 



and the coefficient functions C^{x,^) are defined as : 



(129) 



X — $,-\- is X -\- — is 



(130) 
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In the expression Eq. ()123p for the DVCS amphtude to the twist-3 accuracy, the 
first two terms correspond to the scattering of transversely polarized virtual photons. 
Applying all the definitions to the first two terms of the expression Eq. ()123|1 brings 
us to the twist-2 DVCS amplitude like we had in the previous chapter. The vector 
part 

Ti = \{-gT - jp^){- j^Hp'Mp)Cn + u{p')nu{p)Cu), (131) 

with 

Cn(0 = £^E{x,OC+{x,Odx, (132) 

CuiO = j\{E{x,0 + H{x,0)C+{x,Odx, (133) 
and the axial- vector part 

T2 = --f{-9l' - (^)( - j^u{p'hMp)C2i + u{p')n^Mp)C22), (134) 

with 

C2i(0 = j\E{x,OC-{x,Odx, (135) 

C22(0 = J\ H{x, OC~{x, i)dx. (136) 

This part of the amplitude, containing np and JF^, depends only on the twist-2 
GPDs if, E and if, E and was elaborated in Refs. [45,46]. Comparison with structure 
blocks of Eq. ()1()4|1 gives us the idea of how the structure functions are related to 
the structure functions parametrizing the tensor E{x, ^), vector ll{x^ ^), pseudoscalar 
E{x,^), and axial-vector H{x,^) transitions: 

CuiO YT~^m^^, (137) 
CuiO - CuiO JY^fm, (138) 

C2i(0 ^iT^fpiO, (139) 

Ci2(0 - -^^^fm- (140) 

The generalized parton distributions coincide with the quark distributions at vanish- 
ing momentum transfer and the first moments of the GPDs are related to the elastic 
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form factors of the nucleon through model independent sum rules [1,2] : 

J^ldxH'^{x,U) = Flit), (141) 

J^'^dxE'^ix^U) = F^{t), (142) 

f\xH\x,i,t) = g\{t), (143) 

J^^'dxE'^{x,^,t) = h\{t), (144) 

while the second Mellin moment is closely related to the quark orbital momentum 
contribution to the proton spin [48]. 

The third term in Eq. p23j) corresponds to the contribution of the longitudinal 
polarization of the virtual photon. Defining the polarization vector of the virtual 
photon as 

eiiq) = ^(2^P'^ + ^n^), (145) 



we can easily calculate the DVCS amplitude for longitudinal polarization of the 
virtual photon {L ^ T transition), which is purely of twist-3 : 

(sl)^ H'^'' = ^ £ dx (^J-I C+ix, - ^efAk C-(x, . (146) 

It is therefore seen that this term depends only on new 'transverse' GPDs JF^ 
and jFj^, which can be related to the twist-2 GPDs H, E, H and E with help of 
Wandzura-Wilczek relations given by Eqs. ()124H127j) . Twist-3 part of the amplitude 
could be divided into three separate terms. Two of them 

{u{p')u{p)A±kCsi + u{p')nu{p)A±kC32), 

and 

- ~T (p.q) ^'^^ ' w^y ^^^^^ 

{u{ph5u{p)AnC42 + u{p')n-f5u{p)AnC43), 

are similar to Ti and T2 and give no output to the cross section correction at twist-3 
accuracy. The third one 



(u{p')-fku{p)C33 + ie[u{p')-fa5u{p)C4i 
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with 

C33(0 = f^du!^{E{u,0 + H{u,0) (150) 



^ dx (C+ix, OW+{x, u, - C-{x, OW-{x, u, O) 



C4i(0 = / du^H{u,0 (151) 
dx {C+{x, OW-{x, u, - C-{x, OW+{x, u, O) |, 



-1 
1 



-1 

is the only one that needs to be calculated for the estimate of twist-3 corrections. 
III.3.3 Twist-3 DVCS cross section 

In the leading twist-2 approximation, the amplitude squared falls off as 



Cn - (1 - e)Ci2\' + \eC21 - (1 - e)C22\') ■ (152) 



Next to the leading, twist-3, approximation is suppressed by l/Q"^: 



g4(g2 + t)2(i_^2) 

( - l2e\C33\' + (4eC33 + C41I'). (153) 



Here we used the notation: 

/3^/3, + e = e-4f^. (154) 



.2 • 



Coefficient functions C33 and C41 defined by Eqs. fjl50|) and p5ip are affected by 
the Wandzura-Wilczek transformation. To better understand the properties of the 
WW transformations let us consider two limiting cases of the WW transformation: 
the forward limit ^ — > and the 'meson' limit ^ — 1. In the forward limit we easily 
obtain: 

lim du W+{x,u,0 f(u,n = eix > 0) ^^^^ du 

C-^OJ_i Jx u 

- eix < 0) r ^-^^^ du , (155) 

J -I u 

lim duW_{x,u,^) f{u,C) = 0. (156) 



58 



From the second expression it follows that the W- kernel is a 'genuine non-forward' 
object as it disappears in the forward limit. 

In the limit ^ 1 the generalized parton distributions have properties of meson 
distribution amplitudes [64,65]. In this limit the WW transformations have the form: 

5^1 j-i 2 U-i 1 — u Jx 1 + u J 

In Refs. [50,53,54] it was demonstrated that the twist-3 skewed parton distribu- 
tions in the WW approximation exhibit discontinuities at the points x = Using 
general properties of the WW transformation (jl58p and Eqs. (jl24ll25p one obtains: 

hm /' du{W±{^ + 5, u, - lV±(e - 5, u, = /' ^^^du, 

hm du{W±{-i + 5,u,i) -W^{-i ~ 5,u,0)f{u,i) = (157) 

2 J-i u + ^ 

Here V means an integral in the sense of principal value. We see that for a very wide 
class of functions f{u,^), the discontinuity of the corresponding WW transforms is 
nonzero. This feature of the WW transformation may lead to the violation of the 
factorization for the twist-3 DVCS amplitude. 

Fortunately, some combinations of the distributions JF^ and JF^ are free of discon- 
tinuities due to a certain symmetry of the Eqs. p24|) and ()125|) . For example: 

is free of the discontinuity at x = ^ while its 'dual' combination: 

+ ^^±MP-^p(a;, , (159) 

has no discontinuity at x = — ^. The cancellation of discontinuities in these particular 
combinations of the GPDs ensures the factorization of the twist-3 DVCS amplitude 
on the nucleon. 

One of the non-trivial properties of the generalized parton distributions is the 
polynomiality of their Mellin moments which follows from the Lorentz invariance of 
nucleon matrix elements [1,2,59]: 

f^dxx''H'^{x,0 = K^""^ + hf""^ e + . . . + h'ill e-"', (160) 
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The time reversal invariance requirement [60, 61] leaves only even powers of the 
skewedness parameter ^. This fact implies that the highest power of ^ is + 1 for 
odd (singlet GPDs ) and for even N (nonsinglet GPDs). Furthermore, there is 
a relation between the highest power of ^ for the singlet functions H'^ and E"^ due to 
the fact that the nucleon has spin 1/2 [1,2,61]: 

e^i = -45 . (161) 

The polynomiality conditions p60p strongly restrict the class of functions of two vari- 
ables if^(x, ^) and E'^{x,C,)- In the previous section the polynomiality conditions was 
implemented by using the double distributions [3,14,42]. In this case the generalized 
distributions are obtained as a one-dimensional section (Eq. I103|) of the two- variable 
double distributions F{x,y) (Eq. 11011) [62]. 

It is easy to check that the GPDs obtained by reduction from the double distribu- 
tions satisfy the polynomiality conditions (jl60j) but always lead to h'j^^l = = 0, 
i.e. the highest power of ^ for the singlet GPDs is absent. In other words the 
parametrization of the singlet GPDs in terms of double distributions is not com- 
plete. It can be completed by adding the so-called D-term to Eq. (|103|) [63]: 



jr (x)= dx C S(x + (y-X)F(x,y)dy±9 
Jo Jo 



1-^ 

e 



Note that for both tensor E'^{x,C,) and vector H'^{x,C,) GPDs the absolute value of 
the D-term is the same and has an opposite sign. 

For the quark helicity dependent GPDs H{x,C,) and E{x,^) the D-term is absent. 



III.4 CONCLUSION 

In this chapter, we give the predictions for the DVCS cross section for the kinematics 
reachable by CEBAF at Jefferson Lab. The DVCS amplitude was estimated using 
the model for the non-forward parton distributions from [43]. 

The results for the DVCS part of the total cross section are close to those obtained 
in [57]. 

As was pointed out, the contribution of the interference term to the differential 
cross section is small compared to that of the Bethe-Heitler part alone. This gives the 
idea of how difficult it could be to get the information on real and imaginary parts 
of the DVCS amplitude from the experiment. However, it is a well known fact that 
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the interference term for the cross section of DVCS with the electrons of opposite 
hehcity will give us an estimate of the imaginary part of the DVCS amplitude, while 
the asymmetry of electron-positron cross sections will reveal the real part of the 
amplitude. 

Moreover, at sufficiently large &cm angles the DVCS part of the total cross section 
dominates, but those angles correspond to the large values of \t\. That is why it is 
also necessary to create a "non-zero-f model of the non-forward distributions to 
better estimate the DVCS input from the generalized parton distributions. 
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CHAPTER IV 
CONCLUSION 

In this dissertation two aspects of deeply virtual Compton scattering were studied. 

One aspect considers the small-x DVCS where the energy of the incoming virtual 
photon is very large in comparison to its virtuality. The DVCS in this region can be 
described by the BFKL pomeron. At the same time these kinematics x ~ 10~^ — 
10-^, g2 > 6GeV2, and ~ 1 - SGeV^ could be accessed at HERA [18]. 

It is shown that at large momentum transfer the coupling of the BFKL pomeron 
to the nucleon is equal to the Dirac form factor of the nucleon. It allows us to 
calculate the DVCS cross section and to investigate i-dependence of the amplitude 
in a model-independent way. These approximate calculations of the DVCS cross 
section in the Regge regime in QCD are very timely because all the other predictions 
for the small-x DVCS rely on some model assumptions. 

The algorithm for the calculation of the photon impact factor was elaborated. It 
allows also to estimate the meson production processes amplitude. 

The study showed that the t-dependence of the BFKL pomeron is very important 
and it should be possible to detect at t < 2GeV^. 

Unfortunately, this approximation contains some uncertainties. There is no evi- 
dence that the next ~ (ctslnx)^ term in the BFKL series is as small as we assume. 
Also, we are unable to fix the ag in the leading logarithmic approximation. 

Thus, the first part of this study is devoted to the exploration of the deeply virtual 
Compton scattering process in the Regge regime. It gives an estimation for the cross 
section in soon to be accessible kinematics without any model assumptions. On the 
other hand, it does not reveal any information about the structure of the targeted 
nucleon. 

To study the internal structure of the nucleon, we used the concept of General- 
ized Parton Distributions within the quantum chromo dynamics framework. Those 
distribution functions are accessible through deeply virtual exclusive reactions and 
mainly through the deeply virtual Compton scattering at the kinematical region of 
high virtuality and low momentum transfer t. 

In our study we analyzed the leading twist-2 approximation of DVCS amplitude 
and prepared the algorithm for the estimation of the next-to-leading order, twist-3 
term of the amplitude. In our analysis we neglected the A^^ (twist-3 approximation) 
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FIG. 33: (f) dependence of the beam spin asymmetry A [66]. = 1.25(GeV/c)^, 
xb = 0.19, and -t = 0.19 {GeV/cf. 

which is natural for the parton model. 

DVCS is regarded as the cleanest tool to access the underlying GPDs. Unfortu- 
nately, virtual Compton scattering is always in competition with the dominating BH 
process. 

In our study we created a computer program for precise calculation of the Bethe- 
Heitler process in the required kinematics. 

To exclude the unnecessary input from BH process, the study of asymmetry with 
two electrons of opposite helicities was invoked. It was shown that the main input in 
this case comes from the interference DVCS-BH term, which simultaneously filters 
DVCS observables and magnifies them by the Bethe-Heitler magnitude. 

Numerical algorithms for calculating the DVCS amplitude in twist-2 and twist-3 
approximations were developed. The results for the DVCS cross section at different 
kinematical regions accessible on CEBAF at Jefferson Lab in twist-2 approximation 
were presented. The first measurements of the beam spin asymmetry in the DVCS 
regime were performed at CEBAF at Jefferson Lab [66,71]. As was expected, the 
interference of the DVCS and Bethe-Heitler processes provide a clear asymmetry 
(Eq. ()112|) ). The results of our study based on models from [43] agree in sign, and 
in magnitude to predictions made before [46, 57] and to the experimental results 
represented on Figs. ESI and OH 

It proves that DVCS could be a good tool to access GPDs at relatively low energies 
and momentum transfers. 
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FIG. 34: The new data sample at 5.75 GeV [71] averaged over the entire kinematic 
range of the data. 



With the upcoming CEBAF upgrade to 12 GeV [70] at Jefferson Lab further 
measurements at higher beam energy are expected. Construction of new detectors, 
increase in beam energy, and an improvement of polarization, luminosity, acceptance, 
and precision will broaden the range of the and xb- This will make possible a 
better understanding of nucleon structure. 

Prom the theoretical side of the study, higher order perturbative QCD corrections 
as well as new models of t-dependence of GPDs are expected. 
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APPENDIX A 
BETHE-HEITLER CROSS-SECTION IN DVCS REGIME 

This program was created on the request of experiment aUsts at Jefferson Lab. It 
calculates the cross-section of the Beth-Heitler process at DVCS kinematics. The 
cross-section is calculated in nb/GeV^ with respect to the initial and final electron 
energy and azimuthal and out-of-the-plane angles of the final electron and photon. 

This program uses the Dirac -Ff (i) and Pauli -Ff (0 form factors of the proton as 
a nucleon form factors. 

c Subroutine BH(EO,Ee,ae,ag,f e,fg,cs) calculates the cross 
c section of the e+p->e ' +p ' +gainma reaction, which is dif- 
c ferential with respect to Q**2, Bjorken x, t and out-of- 
c lepton-plane angle f for the process where the photon is 
c emitted from the initial or final lepton ( Bethe-Heitler 
c process) . 

c Input parameters include: 

c EO - initial electron energy, Gev 

c Ee - final electron energy, Gev 

c ae - azimuthal angle of final electron, rad 

c ag - azimuthal angle of emitted photon, rad 

c fe - out-of -plane electron eingle, rad 

c fg ~ out-of -plane photon angle, rad 

c Outcoming parameter is 

c cs - Bethe-Heitler cross section, nb/Gev**4 



c 



Program BetheHeitler 

Pi=3. 141592663 

E0=5 . 75 

Ee=3.75 

ag=20.*Pi/180. 

fe=Pi 

fg=0. 

do 1 i=-24,24 



ae=i*Pi/180. -0.001 

call BH(EO,Ee,ae,ag,f e,fg,cs) 

print 2,ae*180/Pi,cs 

end do 

formate ag=',F10.5,' cs=' ,lpE10.4) 
end 

Subroutine BH (EO , Ee , ae , ag , f e , f g , cs) 
Real M,klq2 
Pi=3. 141592663 
M=.938 

x=EO*Ee* ( 1 . -Cos (ae) ) / ( (EO-Ee) *M) 
y=(EO-Ee)/EO 

Eg=- (-E0* (Ee* ( 1-Cos (ae) ) -M) -Ee*M) / (-E0* (1-Cos (ag) ) -M+ 
+Ee*(l-Cos(ae)*Cos(ag)-Cos(fe-fg)*Sin(ae)*Sin(ag))) 
if (Eg. It. 0) return 
Q2=2*E0*Ee* ( 1-Cos (ae) ) 

t=2* (-E0* (Ee* (1-Cos (ae) ) +Eg* (1-Cos (ag) ) ) + 
+Ee*Eg* (1-Cos (ae) *Cos (ag) -Cos (f e-f g) *Sin(ae) *Sin(ag) ) ) 
klq2=Eg* (1-Cos (ag) ) / (2*Ee* (1-Cos (ae) ) ) 
elch=2.*Sqrt (Pi/137.) 

vph=(x*y**2)/(512. *Pi**4*Q2**2*Sqrt (l+(4*M**2*x**2)/Q2) ) 

bhsq=- ( (2* (2*M**2* ( l+8*klq2**2+t**2/Q2**2+ 
+klq2* (4+4*t/Q2) ) *x**2*y**2+t* (4-2* (2+x+4*klq2*x+t*x/Q2) *y+ 
+ ( (t**2*x**2+2*Q2*t*x* ( l+2*klq2*x) +Q2**2* (2+x**2+8* 
*klq2**2*x**2+4*klq2*x* ( 1+x) ) ) *y**2) /Q2**2) ) * 
*F1 (t) **2) / (x**2*y**2) +4*t* (l+8*klq2**2+t**2/Q2**2+klq2* 
* (4+4*t/Q2) ) *F1 (t ) *F2 (t ) + (t* (M**2* ( l+8*klq2**2+t**2/ 
/Q2**2+klq2* (4+4*t/Q2) ) *x**2*y**2+t* (-2+(2+x+4*klq2*x+ 
+t*x/Q2) *y- ( l+2*klq2*x+t*x/Q2) *y**2) ) *F2 (t ) **2) / 
/ (M**2*x**2*y**2) ) / (klq2*t**2* ( l+2*klq2+t/Q2) ) 

cs=elch**6*bhsq*vph/ (2 . 56942E-6) 

return 

end 

Function Fl(t) 
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f 1=1/ ( 1-t/ . 706) **2* (1-2 . 79*t/ (4* . 938**2) ) / ( 1-t/ (4* . 938**2) ) 
end 

Function F2(t) 

f 2=1/ (1-t/ . 706) **2*1 . 79/ (1-t/ (4* . 938**2) ) 
end 
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